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Summary. — Weak gravitational lensing of background galaxies is a unique, direct 
probe of the distribution of matter in clusters of galaxies. We review several im- 
portant aspects of cluster weak gravitational lensing together with recent advances 
in weak lensing techniques for measuring cluster lensing profiles and constraining 
cluster structure parameters. 
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1. Introduction 

Propagation of light rays from a distant source to the observer is governed by the 
gravitational field of intervening mass fluctuations as well as by the global geometry of the 
universe. The images of background sources hence carry the imprint of the gravitational 
potential of intervening cosmic structures, and their statistical properties can be used to 
test the background cosmological models. 

The deep gravitational potential wells of clusters of galaxies generate weak shape 
distortions of the images of background sources due to differential deflection of light 
rays, resulting in a systematic distortion pattern of background source images around the 
center of massive clusters, known as weak gravitational lensing[l, 2, 3, 4, 5]. In the past 
decade, weak lensing has become a powerful, reliable measure to map the distribution 
of matter in clusters, dominated by invisible dark matter (DM), without requiring any 
assumption about the physical and dynamical state of the system[6, 7]. Recently, cluster 
weak lensing has been used to examine the form of DM density profiles [8, 9, 10, 11, 
12, 13, 14, 15], aiming for an observational test of the equilibrium density profile of DM 
halos and the scaling relation between halo mass and concentration, predicted by ./V-body 
simulations in the standard Lambda Cold Dark Matter (ACDM) modcl[16, 17]. Weak 
lensing techniques have also been used to search for cluster-sized mass concentrations 
projected on the sky [18, 19, 20], allowing us to define samples of "shear-selected" DM 
halos from deep optical surveys[21, 22, 23]. 

In this lecture we briefly review several important aspects of cluster weak gravitational 
lensing. There have been several reviews of relevant subjects: For general treatments 
of gravitational lensing, we refer the reader to Schneider, Ehlers, & Falco[24], Blandford 
& Narayan[25], Refsdal & Surdcj[26], and Narayan & Bartelmann[27]. For a review on 
strong gravitational lensing in clusters, see Hattori, Kneib, & Makino[28]. For a general 
review of weak gravitational lensing, see Bartelmann & Schneider [4]. 

2. Basic Equations of Cosmological Gravitational Lensing 

The local universe appears to be highly inhomogeneous on a wide range of scales from 
stars, galaxies, through galaxy clusters, up to forming superclusters, filaments, and low- 
density voids. The propagation of a light ray is therefore influenced by the gravitational 
field caused by these inhomogeneous mass distributions. The light propagation in an 
arbitrary curved space time is in general a complicated theoretical problem. However, 
most of the astrophysical relevant situations allow for a much simpler description of 
the light ray propagation, which is called the gravitational lensing theory. This section 
reviews briefly the gravitational lensing theory in order to provide a foundation for later 
discussions on cluster weak lensing as well as to summarize the basic equations and 
concepts on cosmological gravitational lensing. 

2T. Gravitational Deflection and Lens Equation. - To begin with, let us consider the 
bending of light rays in the asymptotically flat spacetime caused by a quasi-stationary/spatially- 
isolated mass distribution, assuming that the gravitational field is weak or the deflection 
angle is small. This can be done by solving the null geodesic equation in the Minkowski 
spacetime perturbed with the Newtonian gravitational potential (see, e.g., Ref. [24]). 
The perturbed metric (pi, v = 0,1, .., 3) is then written as 



(1) ds 2 = g^dx^dx" = -(1 + 2^/c 2 )c 2 dt 2 + (1 + 2^>/c 2 )- 1 5 lj dx l dx3 , 




Fig. 1. - Illustration of a gravitational lensing system. The light ray propagates from the source 
S at the position 7] in the source plane to the observer O, passing the position £ in the lens 
plane where it is deflected by an angle a. The angular position of the source S relative to the 
optic axis is denoted by (5, and that of the image I relative to the optical axis is denoted by 
9. The angular diameter distances between the observer and lens, between the lens and source, 
and between the observer and source are Dd , Dds , and D s , respectively. 



where * is the Newtonian gravitational potential, 5ij — 1,2,3) is the Kronecker 
delta, and c is the speed of light. We consider the metric given by equation (1) to be 
the sum of a background metric g$ and a small perturbation h pv : g pv = g$ + h pv 
with \h pu \ -C 1. Expanding to the first order in <F/c 2 we break up the metric in equation 
(1) into a sum of the background Minkowski metric g$ = r)^ = diag(— 1, 1, 1, 1) and a 
perturbation of the form: h pv = diag(-2>F/c 2 , -2*/c 2 , -2*/c 2 , -2*/c 2 ). As usual, g^ v 

and g {b ^ u are defined by g^g pv = 5% and g (b)tip gp b J = 5£. Then, to the first order of h, 
g» v = g^ v - h» v , where W v is defined by W v = g^Pg^ va h p(T . 

The propagation of light rays is described by the null geodesic equation: 



(2) 
(3) 

(4) 



k ~ dX ' 
= g^i/k^k , 

dkf * - -F^ 1 h v h x 
~~d~X ~ Kxk k ' 



where fc M is the 4-momentum of the light ray, A is the affine parameter, and T^ x de- 
notes the Christoffel symbol: T^ p = (l/2)g AlA {g\ VtP + g\ p , v — 9u P ,x)- In the background 
Minkowski spacetime with gffl = r/ pi/ , T^p^ = 0. For a light ray moving into the pos- 
itive x-direction in the background spacetime, the photon 4 momentum k^P and the 
corresponding unperturbed orbit x^^ = (ct,x) are simply given as k^^ — dx^^/dX = 
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(1,1,0,0) and = (A, A, 0,0). 

Now we consider the light ray propagation in the perturbed space time. Without loss 
of generality, we can take the deflection angle to lie in the xz plane (i.e., y = x 2 = 0). 
In the weak field limit (I'J'/c 2 ! <C 1) the impact parameter b of the incoming light ray is 
much greater than the Schwarzschild radius of the deflector, b ^> 2GM/c 2 . The linearized 
null geodesic equation takes the form: 



(5) 
(6) 

(7) 



jfe"(A) = kS b ^{X) + 5k fl (X), 

= h^k^k^ + 2g$k^5k u , 
d(<5fc") 



dX 



-2T^kV>>6k x - 6I* x k<V v kW\ 



Note that the first term on the right-hand side of equation (7) vanishes since 
The perturbed Christoffel symbol is 5T£ p = {l/2)rf x (hx„, p + h Xp ,„ - K pA ) + 0(h 2 ). 
Choosing the boundary condition <5/c A '(— oo) = 0, we integrate the linearized geodesic 
equation (7) to obtain the 4 momentum of the light ray in the out state as: 



(8) 



5jfe"(+oo) = < 



2f^d\d ]] */c 2 ( M = 0) 
(M=l) 

2f+™dXd ± */c 2 ( M = 2) 



In the small angle scattering limit 
obtained as 



|*/c 2 | <C 1), the bending angle a of the light ray is 



(9) 



a : 



fc ± (+oo) 
fcll(-oo) 



2 f + °° 
~# 1 ' 

L J — oo 



dXd±^(X). 



For (j, = 1 (the x" -component), <5/c"(A) — const. = 0, so that fell (A) = 1 and xll = 
A + 0(\h\). This suggests a simple approximation for the bending angle consistent with 
the weak field limit: Using dX — dx^ + 0(\h\) and integrating along the unperturbed 
path, we have 



(10) 



a 



9 r+oo 

— 2 / dx 11 V ± ^(x ll ,x±). 

c J — oo 



This is known as the Born approximation. This yields an explicit expression for the 
bending angle as a w AGM/(bc 2 ) w 1.75"(M/M )(6/ J R )~ 1 . General relativity gives 
a deflection angle twice as large as that Newtonian physics would provide. Eddington 
confirmed the prediction of general relativity from the measurement of the deflection 
of the starlight during a solar eclipse. This solar deflection angle is verified within ~ 
0.1%[29]. 

Finally, we find from the null condition (equation [3]) that Sk°(\) = — 2\E , (A) + 0(h 2 ), 
or dt/dX = 1 — 2*(A) + 0(h 2 ) > 1. The gravitational time delay Ai grav , with respect to 
the unperturbed light propagation, is thus given by 



(11) 



cAt 



grav 



/+oo 
-oo 



dA*(A)/c 2 
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Let us consider a situation as illustrated in Figure 1: A light ray propagates from a 
source S at the position rj in the source plane to an observer O, passing the position £ in 
the lens plane where it is deflected by a bending angle a. Here the source and lens planes 
are defined as planes perpendicular to the optical axis at the distance of the source and 
the lens, respectively. The exact definition of the optical axis does not matter since the 
angular scales involved in a typical lensing situation are very small. The angle between 
the optical axis and the true source position is (3, and the angle between the optical axis 
and the image / is 9. The angular diameter distances from the observer to the deflector, 
from the observer to the source, and from the deflector to the source, are denoted by 
Dd, D s , and Dd s , respectively. From Figure 1, we find the following geometrical relation: 
t] = (D s /Dd)£ — Dd s 6t(£)- Equivalently, this is transformed to the relation between the 
angular source and image positions, (3 — rj/D s and 9 = £/Dd'- 

(12) (3 = e+^dt{D d 6) =9 + a{6), 

where we have defined the reduced deflection angle a(9) in the last step. Equation 

(12) is called the lens equation, or ray-tracing equation. The lens equation is in general 
non- linear (with respect to 9) , so that it may have several solutions 9 for a given source 
position /3, corresponding to the multiple imaging of a single source on the sky. 

2'2. Cosmological Lens Equation. - Here we turn to the cosmological lens equation 
that describes the light propagation in an expanding, locally-inhomogeneous universe. 
There are various approaches to derive the cosmological lens equation[30, 31, 32, 33]. 
Here we follow the approach by Futamase [33] based on the linearized null geodesic 
equation as introduced in §2T. 

Consider the following perturbation of the Friedman-Lcmaitre- Robertson- Walker (FLRW) 
metric: 

(13) ds 2 = a 2 {fi)d~s 2 = a 2 (r ] )g^dx f *dx l/ , 

(14) = a 2 [-(1 + 2^)d V 2 + (1 - 2*) {d X 2 + r 2 { X ){d6 2 + sin 2 0d(f> 2 )}] , 

where (x M ) = (77, x, 9, <f>), 77 = c f'dt' /a(t') is the conformal time, and r(x) is the comoving 
angular diameter distance defined by 

A"- 1 / 2 sin (x/^/k) K > 

The curvature parameter K is expressed with the present-day total density parameter of 
the Universe Oo as K — (Oo — ^-)Hq/c 2 . Since the structure of a light cone is invariant 
under the conformal transformation, we work in the conformally-related spacetime with a 
null geodesic given by ds 2 . The metric can be rewritten in the form of = gffl+h^, 
as a sum of the background FLRW metric and a small perturbation (\h\ <C 1). The 
cosmological Poisson equation relates the Newtonian gravitational potential $ to the 
density perturbation field, Sp m , as: 



(15) r( X ) = { 



(16) V 2 *( X ) = AnGa 2 5 Pm = ' 



m u m 
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where 5 m — Sp m /p is the overdensity with respect to the mean background density of 
the universe, fi m is the present-day matter density parameter, and p is the mean cosmic 
matter density, p = fl m a~ 3 /0 C rit,o — a _3 (3i?o^m)/(87rG). The physical implication of 
equation (16) is that the amplitude of >J/ is related to the amplitude of S as |^/c 2 | ~ 
(3Cl m /2)(l I 'Lh) 2 '(d m / a) where I and Lh = c/H denote the characteristic comoving scale 
of the density fluctuation and the Hubble radius, respectively. Therefore, assuming the 
standard power spectrum of the density perturbation, we can safely conclude that the 
degree of metric perturbation is always much smaller than unity, i.e., |*/c 2 | <C 1, even for 
highly non-linear perturbations with |i5 m | ^> 1 on small scales of I <C Lh(— 3Gpcft. _1 ). 

Following the prescription given in §2T, we solve the null geodesic equation (equations 
[2], [3], [4]) on the perturbed FLRW metric (13). Choosing the origin of the coordinate 
system at the observer's point and backward ray-tracing from the observer A = to the 
source A = A s along the photon path, we obtain the 4-momcntum of the light ray in 
the background FLRW spacetime as = (—1, 1, 0, 0), and the corresponding orbit as 
x (b ^(X) = (-A,A,0/,0 7 ) with 6»/ < 1, where (0j, 0j) = (0(0), 0(0)) denotes the angular 
direction of the image on the celestial sphere. The comoving angular diameter distance 
r(x) defined by equation (15) can be parametrized by the affine parameter along the 
photon path as r(x(X)) = r(A). The null geodesic equation for the perturbative part is 
given by equations (5), (6), (7), and can be formally solved as 

(17) 5k"(X) = ~^-r [^dX'd^W/c 2 (p = 9, 0), 

r (A) Jo 

where = (^,e, sin -2 9^^). Inserting this in equation (5) and integrating the angular 
part yields the following equation: 

(18) , Ss ^. ) = 9f -l|\ A l(^)^ W , 

where quantities with the subscript S denote those defined for the source. Now we 
consider a small patch of the sky over which the curvature of the sky is negligible (flat-sky 
approximation). Then, one can define locally flat-sky Cartesian coordinates 9 = (9,9(j>) 
around the line-of-sight of interest (9 <C 1). Defining 13 = 6s and = Oj, equation (18) 
is written in this coordinate system as 

/•Source 

(19) 13-0= da = a ( X s), 

J Observer 

(20) a( Xs ) = ~ [ XS d\ T{Xs ~ X) V ± nx(X)); x(X) - x( b \X) + 5x(X), 

where Vj^ is the (comoving) transverse derivative, Vj_ = r~ 1 (X)(de,9~ 1 dj > ), and the 
integral is performed along the perturbed trajectory x(X) = x^(X) + Sx(X) with X = x + 
0(| ^/c 2 1). This is a general expression of the cosmological lens equation: as long as the 
approximations adopted are valid (see Refs. [33] and [34] for more details), equation (19) 
can be applied to various lensing phenomena, including the multiple deflections of light 
rays, strong and weak lensing by clusters of galaxies, and the weak lensing by the cosmic 
large-scale structure continuously distributed between the source and the observer. Note 
that the cosmological lens equation (19) is derived with the standard angular diameter 
distance in FLRW spacetime without employing the thin lens approximation. 
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2'3. Cluster Gravitational Lensing. - Now, let us turn to the case of gravitational 
lensing by clusters of galaxies. Clusters are the largest gravitationally-bound systems 
observed in the Universe, with typical masses of M ~ 10 14 - 10 15 M Q . In the context 
of the standard structure formation scenario, clusters, as the most massive collapsed 
structures, correspond to rare peaks in the primordial density field. Clusters produce 
deep gravitational potential wells, and act as efficient gravitational lenses. The large- 
scale structure in the Universe also affects the propagation of light rays from distant 
sources, producing small but continuous transverse excursions of light rays along the 
light path. In the study of cluster gravitational lensing, it is often assumed that the total 
deflection angle a is dominated by the gravitational potential of the cluster itself, and 
that contributions from the cosmic large-scale structure and multiple deflections by other 
non-linear objects are negligible. Assuming that the light propagation is approximated 
by a single scattering event by the cluster (single lens approximation) and that a light 
deflection occurs within a sufficiently small region (xd — Ax/2, Xd + Ax/2) compared to 
the relevant angular diameter distances (thin- lens approximation) , we have the deflection 
angle by a single cluster as 



where D s — a{Xs)f(Xs) and Dd. s = a(Xs) r (Xs — Xd) are the angular diameter distances 
from the observer to the source, and from the deflector to the source, respectively, and 
r{xd)& is the comoving transverse vector on the lens plane. In a cosmological situation, 
the angular diameter distances Dij between the planes i and j (zi < Zj ) are of the order 
of the Hubble radius, Lh = c/Hq w 3Gpc/i _1 , while physical extents of clusters are 
about 2R = 2 — 5Mpc/i _1 . Therefore, one can safely adopt the thin-lens approximation 
in cluster gravitational lensing. 

In actual observations, large scale structure along the line-of-sight also contributes to 
the lensing signal, and consequently affects the measurements of cluster mass properties. 
The level of uncertainties on cluster lensing measurements due to large scale structure 
can be assessed by numerical simulations (see Ref. [35]). For a massive cold dark matter 
(CDM) halo with M 2 oo = 10 15 M Q /i _1 , uncertainties in the mass measurement could 
reach 16 — 18% at intermediate redshifts of = 0.2 — 0.3[35]. See also Appendix A 
for the multiple lens equation based on a discretized version of the cosmological lens 
equation. 

Then, we introduce the effective lensing potential ip(8) defined as 



where Dd is the angular diameter distance from the observer to the deflector, Dd = 
a(xd) r (Xd)- In terms of ip(0), the lens equation is expressed as 



where Ve = rVi = {d$,6 1 d r j > ). Table 1 gives examples of analytic axially-symmetric 
lens models based on spherically-symmetric mass distributions, such as the point mass, 



(21) 




(22) 




(23) 
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Table I. - Examples of analytic lens models based on the spherical symmetric mass distribution. 
The 3D density profile p(r), the lensing convergence k(6), and the mean convergence R{6) inside 
the angular radius 9 are given. Here 'T' and 'R' indicates that the lens has a tangential and a 
radial critical curve, respectively. 



Model 



p(r) 



K(0) 



Remarks 



Point mass 



MS'Uf) 



A-kGM D d D dB a * 

— d— 5d{6) 



AGM DdDde 1 
c 2 D a e 2 



SIS« 



2-wGr 2 



\ c J D s 6 V c J D s 6 



T 



isc< 2 > 



P(0) 



1 + (r/r c ) 2 



k(0) 



V 1 + ( W 2 



2«(0) 



T, R 



NFW (3) 



pcrit^c 



(r/r s )(l + r/r s ) 2 



Ksf(o/e a ) 



2n s g(e/e s )/(e/e s y 



T, R 



Notes. 

(1) Singular isothermal sphere: o v represents the isothermal ID velocity dispersion. 

(2) Isothermal sphere with a finite core radius: r c represents the core radius, and 6 C is its angular 
radius, 6 C := r c /Dd\ k(0) := £ m (0)/£ cr it with E m (0) = 7rp(0)r c . The ISC produces two critical 
curves if, and only if, k(0) > 1. 

(3) Navarro, Frenk, and White (NFW[37]) universal profile of cold dark matter (CDM) halos: 5 C , 
Pcrit, and r s represent the characteristic overdensity of the CDM halo, the critical density of the 
Universe, and the scale radius, respectively. k s and 8 S are then defined by k s := 2<5 c pcritr s E ( 
and 9 B := r s /Dd, respectively. The functions f(x) and g(x) are defined as follows [38] [39] 



i 

crit 



(24) 



(25) 



' i^(-i + 7= arct W^) (x<1) 



= in (|) 



+ 1- 



■ arctan 



+ < 



= arctanh 



arctan , , 
I v x+1 




the singular isothermal sphere, and the NFW models. For a more complete catalog of 
mass models for gravitational lensing, see Keeton [36] . 

2'4. Lensing Jacobian Matrix. - The local properties of the lens mapping are described 
by the Jacobian matrix A : 



(26) 
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This symmetric 2x2 Jacobian matrix can be decomposed into the following form: 

(27) A a f} = 5 a f3 ~ 1p,a0 = (1 - K)5 a p - 7i<7 3 - 72^1 

where k(0) is the lensing convergence field defined by 



(28) 



« = \ (V'.u + ^,22) = ^A ^(6»), 



responsible for the trace-part of the Jacobian matrix, "y a (6) (a = 1, 2) are the components 
of the complex shear field j(0) :— 71 (0) + i-f 2 (0), defined as 



(29) 



71 = \ 0,11 - i>,1l) ; 72 = ^ ("042 + 0,21) = 0,12, 



and cr a (a = 1,2,3) are the Pauli matrices that satisfy <J a Vb — i^abc^c- Equation (28) 
can be regarded as the two-dimensional Poisson equation. Then, the Green's function 
(or the propagator) in the infinite domain (1Z 2 ) is A^ 1 (6, 0') = In \6 — 0'\/(2tt), so that 
ip{0) = (1/tt) f n2 d 2 9' \n(0 — 0')n{0'). The explicit representation of the lens Jacobian 
matrix is 



(30) A(0) 



1 - K - 71 -72 
-72 1 - K + 71 



(1-K) 



1 
1 



71 72 

72 -71 



and it has two eigenvalues A± = 1 — k ± |7|. In Figure 2 we illustrate the effects of 
the lensing convergence k and the gravitational shear 7 on the shape and size of an 
infcnitesimal circular source. 

Convergence alone 



Gravitational Lensing 




Source 




Convergence + Shear 



Fig. 2. - Illustration of the effects of the convergence k and the shear 7 on the shape and size of 
a hypothetical circular source. The convergence acting alone causes an isotropic magnification 
of the image (dashed circle), while the shear deforms it to an ellipse. 
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2'5. Lensing Convergence. - The lensing convergence k is essentially a distance weighted 
mass overdensity projected along the line-of-sight. We express n(9) due to cluster grav- 
itational lensing as 

(3D m = 1 4, (,„ . « (^^_) - K ^ffi, „ _ 

where S m = / c£Z (p m — p) w / p m is the surface mass density field of the lensing cluster 
projected on the sky, and S cr i t is the critical surface mass density of gravitational lensing, 

where dij = D^/Lh is the dimensionless angular diameter distance between the planes 
i and j; £ C rit depends on the lensing geometry (zd, z s ) and the geometry of the Universe, 
e.g., (n m ,fl A ,H ). 

For a given lens redshift Zd, the geometric efficiency of gravitational lensing is de- 
termined by the distance ratio, Dd s /D s , given as a function of the background redshift 
z s and the cosmological parameters. Figure 3 compares Dd s /D s as a function of z s 
for various sets of the lens redshift and the cosmological model. In order to convert 
the observed lensing signal into physical mass units, one needs to evaluate the depth 
of background sources (i.e., the source redshift distribution, or its moments). Figure 4 
shows S~| t oc Dds/D s as a function of lens redshift Zd for three different source redshifts, 
z s = 0.8, 1.0, 1.2. In typical optical imaging observations down to a magnitude limit of 
R c <~ 26 ABmag, the median depth of background galaxies is about z s <~ 1. When the 
lens redshift is Zd £ 0.2, £ cr ;t depends weakly on the source redshift, so that a precise 
knowledge of the redshift distribution of background galaxies is not crucial (see for de- 
tails, e.g, Refs. [40] and [13]). On the other hand, this distance dependence of the lensing 
effects can be used to constrain the cosmological redshift-distance relation by examining 
the geometric scaling of the lensing signal as a function of the background redshift (see 
Refs. [41] and [42]). 

2'6. Quadrupole Shape Distortion: Gravitational Shear. - The deformation of the 
image for an infinitcsimally small background source (d/3 — > 0) can be described by the 
inverse Jacobian matrix A~p = (A~ 1 ) a f3 of the lens equation (a, (3 — 1,2). In the weak 
lensing limit, 

(33) A-^(l + K)S al3 + T al3 , 
where 

(34) r Q/3 = (d a d [3 - 5 af s^Ag^ ip(0) = 0-371 + CT172 

is the symmetric, trace-free 2x2 shear tensor [4] [43]. The first term in equation (33) 
describes the isotropic light focusing (area distortion) in the weak lensing limit, while 
the second term induces an asymmetry in lens mapping; the shear is hence responsible 
for the shape distortion. Note that both the convergence and the shear contribute to the 
shape/area distortions in general (non-weak) cases. 
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Fig. 3. - Distance ratio Dd s /D 3 as a function of the source redshift z 3 for various sets of the 
lens redshift Zd and the cosmological parameters (f2 m ,f2A). Dds/D s is plotted for three lens 
rcdshifts z d = 0.2, 0.5, 0.8 (from left to right), and for three sets of the cosmological parameters: 
Q A ) = (1, 0), (0.3, 0), and (0.3, 0.7). 



2"7. Area Distortion: Gravitational Magnification. - The determinant of the Jacobian 
matrix is given as det*4 = (1 — k) 2 — \j\ 2 - In the weak lensing limit where \k\, 1 71 1, 1 72 1 <C 1, 
det*4 « 1 — 2k. Gravitational lensing describes the light ray deflection in the weak 
field limit (I'J'/c 2 ! <C 1). The surface brightness of a background source is unchanged 
under gravitational lensing (Liouville's theorem). The flux magnification in gravitational 
lensing is due to the light-ray focusing that causes the area distortion: 5Q 1 = fi5Q s . The 
magnification is hence given by taking the ratio between the lensed to the unlensed image 
solid angle as fi = SQ 1 /Sfl s — 1/detA: 

1 _ 1 

(35) M ~dct~4 ~ (1_ K )2_| 7 |2- 

In the weak lensing limit, the magnification to the first order is 

(36) n rj 1 + 2k. 

Thus, for an image at k(6) = 0.1, the corresponding magnitude change is Am w 
-(5/2) log 10 ( M )~ -0.20. 
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Fig. 4. - Geometric scaling of the weak lensing signal. The lower panel shows the inverse 
of the critical surface mass density of gravitational lensing, £~ it (z<i, z s ), as a function of lens 
redshift Zd for three different source redshifts, z s = 0.8, 1.0, 1.2 {dashed, solid, and dotted- dashed, 
respectively), demonstrating the geometric scaling of gravitational lensing signal. The top panel 
shows the relative lensing strength E~; t (zd, z s )/S~J t (zd, z s = 1.0) as a function of lens redshift 
Zd normalized with respect to the source at z s = 1. For lensing clusters at low redshifts Zd, ^Jcrit 
depends very weakly on the background redshift z 3 , so that the uncertainty in z 3 of background 
galaxies is less important in the lensing-based cluster mass determination. This figure is taken 
from Ref. [7]. 



3. Weak Gravitational Lensing 

In this section, we assume for simplicity that the lensing fields are subcritical every- 
where, i.e., detA(6) > 0. 

3'1. Weak Lensing Mass Reconstruction. - For convenience we define the complex 
gradient operator d :— d\ + idi that transforms as a vector, & = de %v , with ip being the 
angle of rotation relative to the original basis[44, 45, 46]. In terms of the effective lensing 
potential ip(0), the lensing convergence is given as 

(37) k(6) - \d*di,{0), 

where * denotes the complex conjugate, and dd* = Vg is a scalar or a spin-0 operator. 
Similarly, the complex shear 7 = 71 + 272 = |7|e 2 ^ is given as 



(38) 
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Offset RA (arcmin) 

Fig. 5. - The projected mass distribution k(0) of A1689 reconstructed using the weak gravi- 
tational shear field j(6) measured from a color/magnitude-selected sample of red background 
galaxies registered in deep Subaru imaging observations. Overlaid up on the image is the recon- 
structed spin-2 gravitational shear field 7(6). A stick with the length of 10% shear is indicated 
in the top right corner. The north is to the top, and the east is to the left. This figure is based 
on the Subaru weak lensing data presented in Ref. [9]. 



where T>q = dd/2 = (d 2 — 9|)/2 + idid2 is a spin-2 operator, which transforms as 
T)' g = 2Ve 2lv under a rotation of the basis axis by (p. 

Since the 7- and K-fields are linear combinations of the second derivatives of ?/>(#), 
7(6) and k(6) are related with each other by (see Refs. [47] and [43]) ( x ) 

(39) Agn(6) = d*d*j = 2V* gl (0) 

Thus, the shear-to-mass inversion formula can be formally obtained as 

(40) k{6) = A," j, [5*5*7(0')] = 2^A^ [7(C)] • 

For the two-dimensional Poisson equation, the Green's function (propagator) A _1 (<?, 0') 
in the infinite domain (1Z 2 ) is A _1 (0, 9') = In \9 — 9'\/ (2ir), so that equation (40) can be 
solved to yield the following non-local relation between k and 7 (see Ref. [2]): 

(41) k(0) = ~ [ d 2 9' D*{9-9') 1 {9'), 



( 1 ) An alternative but equivalent expression is AgK,(0) — d a d lS V a is(0). 




Fig. 6. - Comparison of the surface mass density field and the cluster galaxy distributions 
in C10024+1654. Left: Dimensionless surface mass density field, or the lensing convergence 
k(0) = E m (0)/E cr it, reconstructed from Subaru distortion data. Middle: Observed surface 
number density distribution E n (#) of _Bi? c z'-selected green galaxies, representing unlensed clus- 
ter member galaxies. Right: Observed i? c -band surface luminosity density distribution E;(#) of 
the same cluster membership. The solid circle in each panel indicates the cluster virial radius of 
r v i r — 1.8Mpc/i _1 at the cluster redshift of z — 0.395. All images are smoothed with a circular 
Gaussian of FWHM 1.4arcmin. Also overlaid on the E n (0) and £;(#) maps are the k(8) field 
shown in the left panel, given in units of 2a reconstruction error from the lowest contour level of 
3<r. The field size is 26' x 22'. North is to the top, east to the left. Figure taken from Ref. [40]. 



where D(0) is the complex kernel defined as 



(42, iw-^-w. g-y** -^,. 

Similarly, the complex shear field can be expressed in terms of the lensing convergence 
as 

(43) 7 (0) = I/ d 2 9' D{6-e')K{6'). 

K J-R? 

In a practical application, it is computationally fast to work in Fourier domain [48], by 
using the fast Fourier transform (FFT). By taking the Fourier transform of equation 
(39), we have a mass inversion relation in Fourier space as 

(44) m ^'^'T^ m (MO) 



where k(k) and 7 (ft) are the Fourier transform of the k(6) and j(9) — 7l (0) + 172(0) 
fields, respectively. 

Figure 5 shows the two-dimensional mass distribution in the central 22' x 22' region 
of A1689 at Zd — 0.183 reconstructed from the weak shear field[9] , measured from a 
sample of blue+red background galaxies registered in deep Vi' images taken with the 
Suprime-Cam[49] on the Subaru telescope. Also overlaid up on the image is the gravi- 
tational shear field j(6) of the red background galaxies, revealing a coherent tangential 
pattern around the cluster center. In the left panel of Figure 6 we show the mass map for 
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CL0024+1654 (zd — 0.395) reconstructed from Subaru distortion data of -Bi? c -2'-selected 
blue+red background galaxies. Also compared in Figure 6 are member galaxy distri- 
butions and in the cluster, Gaussian smoothed to the same resolution of 
FWHM = 1.41arcmin. Overall, mass and light are similarly distributed in the cluster. 

Adding a constant mass sheet to n in equation (43) does not change the shear field 
7(0) which is observable in the weak lcnsing limit, leading to the so-called mass-sheet 
degeneracy based solely on shape-distortion measurements [4, 5]. As we shall see in §3'2, 
the observable quantity is not the gravitational shear 7 but the reduced shear, 

(45) 9 (0)= ^ 



k{6) 

in the subcritical regime where det.4 > (or 1/g* in the negative parity region with 
detA < 0). We see that the g-field is invariant under the following global transformation: 

(46) k(9) -> \k(6) + 1 - A, 7 (6») -> A 7 (0) 

with an arbitrary scalar constant A 7^ 0[3]. This transformation is equivalent to scaling 
the Jacobian matrix A(6) with A, A(6) — > \A(9). This mass-sheet degeneracy can be 
unambiguously broken by measuring the magnification effects, because the magnification 
fi transforms under the invariance transformation (46) as 

(47) M (0) -+ AVW- 




x (pixel) x (pixel) 

Fig. 7. - Example of the anisotropy point-spread function (PSF) correction (Subaru i' data of 
A1689). The quadrupole PSF anisotropy field as measured from stellar ellipticities before and 
after the PSF anisotropy correction. The left panel shows the raw ellipticity field of stellar ob- 
jects, and the right panel shows the residual ellipticity field after the PSF anisotropy correction. 
The orientation of the sticks indicates the position angle of the major axis of stellar ellipticity, 
whereas the length is proportional to the modulus of stellar ellipticity. A stick with the length 
of 5% ellipticity is indicated in the top right of the right panel. This figure is taken from Ref. [9] . 



3 '2. Weak Lensing Distortion Observables. - In a moment-based approach of weak 
lensing shape measurements (see Kaiser, Squires, & Broadhurst[50]; hereafter KSB) we 
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Fig. 8. - Stellar ellipticity distributions before and after the PSF anisotropy correction, corre- 
sponding to Figure 7. The left panel shows the raw ellipticity components (e^e^) of stellar 
objects, and the right panel shows the residual ellipticity components (8e*,5e2) after the PSF 
anisotropy correction. This figure is taken from Ref. [9]. 



use quadrupole moments Q a p, (a, /3 = 1,2) of the surface brightness distribution 1(9) of 
background images for quantifying the shape of the images: 

(48) QaP = J<P0qj[I(0)] ' 

where qj[I(0)] denotes the weight function used in noisy shape measurements and A9 a = 
a — 9a is the offset vector from the image centroid.( 2 ) The complex ellipticity e is then 
defined as 

(49) e = Oil ~ Q 22 + 2i Q 12 

Qn + Q22 

The spin-2 ellipticity e transforms under the lens mapping as 

(50) e(*) = e ~ 2 9 + 9 2 e* 

{ 1 1 + \g\ 2 - 2Re [ge*] ' 

where quantities with subscript "(s)" represent those of (unlensed) intrinsic background 
sources, and g = 7/(1 — k) is the spin-2 reduced shear. Since e is a non-zero spin quantity 
with a direction dependence, the expectation value of the intrinsic source ellipticity e' s ^ is 
assumed to vanish: (e^) = 0. Schneider & Seitz[3] showed that = (x^) is equivalent 



( 2 ) A practical implementation of the KSB method is achieved by the IMC AT package developed 
by Nick Kaiser. Note that IMCAT measures the shape moments with respect to the peak position 
rather than the centroid. 
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to 



(51) = 2^ w ii 



K(J fl e;) 



where <5 ff is the spin-2 complex distortion 5 g = 2g/(l + |g| 2 )[3], ej is the complex ellipticity 
for the jth object, and Wj is a statistical weight for the jth object. Thus, for a intrinsically 
circular source with e< s ) = 0, we have 

(52 > 

In the weak lensing limit (|k|, | — ^ | <C 1), equation (50) reduces to w e — 2g w e — 
27. Assuming the random orientation of the background sources, we average observed 
ellipticities over a sufficient number of images to obtain 

(53) (e) « 2g « 2 7 . 

For g( tru °) = 0.1, the weak lensing approximation (53) gives </ cst ) = 0.099, or a negative 
bias of —1%. For g( true ) = 0.4 in the non-linear regime, equation (53) gives g( cst ) = 0.345, 
corresponding to a negative bias of —14%. 

For a practical application of weak lensing shape measurements, we must take into ac- 
count various observational effects such as noise in the shape measurement due to readout 
and/or sky background and the dilution of the lensing signal due to the isotropic/anisotropic 
point-spread function (PSF) effects. Thus, one cannot simply use equation (53) to 
measure the gravitational shear field. KSB took into account explicitly the Gaussian 
weight function in calculations of noisy shape moments and the effect of quadrupole 
PSF anisotropy, as well as isotropic PSF smearing, and derived in the limit of linear 
anisotropies relevant transformation equations between unlensed (intrinsic) and lenscd 
(observed) ellipticities. In the limit of linear anisotropies, the transformation equation 
between the intrinsic and the observed complex ellipticities is formally expressed as 

(54) e a = e^ + (C) ap gp + (C^) ap qp, 

where q a is the spin-2 PSF anisotropy kernel, and Cs are linear response coefficients for 
the spin-2 anisotropy fields (g a and q a ), which can be calculated from observable weighted 
shape moments of galaxies and stellar objects[4, 50, 51]. In practical observations, the 
spin-2 PSF anisotropy q{6) can be measured from image ellipticities of foreground stars, 
for which both of e^ s ^ and g vanish: qp{9) = {C q )~^e*^. As an example, we show in 
Figure 7 the quadrupole PSF anisotropy field as measured from stellar ellipticities before 
and after the anisotropic PSF correction using the Subaru i' data of the cluster A1689[9]. 
Figure 8 shows the distribution of stellar ellipticity components before and after the PSF 
anisotropy correction, corresponding to Figure 7. 

Assuming that the expectation value of the intrinsic source ellipticity vanishes, we 
have the linear relation between the averaged image ellipticity and the reduced shear as 

(55) g a *((C°)-l(e-C«q)e). 
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In the context of the KSB formalism, the linear response C a a p to the gravitational shear 
is often denoted as P a p (or P^p)- A careful calibration of P 9 is crucial for accurate mea- 
surements of the weak lensing signal. See Refs. [52] and [53] for more detailed discussions 
on the shear calibration issues. 

3'3. E/B Decomposition. - In general, the shear tensor field T a p(6) = 71(0)03 + 
72(^)01 that describes the spin-2 quadrupole shape distortions can be expressed as a 
sum of two terms, corresponding to two degrees of freedom, by introducing two scalar 
functions <&e(0) and $b(0) (see Ref. [43]) as 

(56) v a p{9) = ( £ ) =r { a E p\e) + r%\e), 

(57) rg } = (d a 8p - Sap^Ae^j * E {0); T { a B J = 1 (e^d a d 7 + e^dpdj $ B (0), 

where e a p is the 2x2 antisymmetric tensor, en = 622 = 0, £12 = —£21 = 1- The first 
term associated with <&e is a gradient or scalar E component, and the second term with 
$b is a curl or pseudoscalar B component. The 71 and 72 are then written in terms of 
<f>£ and $ B as 

(58) 71 = +rn = -r 22 = ^ ($e,h - $£,22) - $s,i2 

(59) 72 = r 12 = r 21 = $ B)12 + i ($ B)11 - $ Bi 22) • 

As we have seen in §3'1, the spin-2 7 field is coordinate dependent, which transforms as 
7' = -fe 2l{p under a rotation of the basis axis by ip. The E and B parts can be extracted 
from the shear tensor by applying the Vg operator: 

(60) 2V%k e = Vfo E = 2d a d?r a p, 2V|k b ee = 2e Q ^ Q 9 7 r^ 7 , 

where we have defined the E and S fields, ke = (l/2)Ae$B and k b = (l/2)Ae$s, 
respectively. This technique is called as the S-B mode decomposition. The above equa- 
tions tell us that the relations between the S/B-fields and the spin-2 field are intrinsically 
non-local. 

Remembering that in the case of gravitational lensing the shear tensor is given as 
T a p = (d a dp-8 a pA e /2)il>(0), we can identify $ E (0) = ip{0) and $ B (0) = 0. Hence, for 
a lensing-induced shear field, the -E-mode signal is related with the convergence n, i.e., the 
surface mass density of the lens, while the B-mode signal is identically zero. We note that 
gravitational lensing can give rise to -B-modes, for example, when multiple deflections 
of light rays along the light path are involved. However, these B modes arise at higher 
orders and the B-mode contributions coming from multiple deflections are suppressed by 
a large factor compared to the _E-mode contributions. In practical observations, intrinsic 
cllipticities of background galaxies also contribute to the gravitational shear estimate, 
7 cst . Assuming that intrinsic ellipticities have random orientations in projection space, 
such uncorrelated cllipticities will yield statistically identical contributions to the E- and 
B modes. Thus the B-mode signals serve as a useful null check for systematic effects 
(e.g., residual PSF anisotropies). 
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Now we turn to the issue of E/B-mode reconstructions from the spin-2 shear field. 
Rewriting equations (60) and (60) in terms of the complex shear 7, we find 

(61) A s k e = 5ft (2^7) , 

(62) A gKB = 3 (22^7) = -R (2^7) . 

Defining k = ke + in E , we see that equations (61) and (62) are identical to the mass 
inversion equation (39). Therefore, the -B-mode convergence Kb can be simply obtained 
as the imaginary part of equation (41), which is expected to vanish for a purely weak 
lensing signal. Furthermore, the second equality of equation (62) indicates that the trans- 
formation 7'(0) = ij(d) ("/[ = —72,72 = 7i) is equivalent to an interchange operation of 
the E and B modes of the original maps: n' E {0) = —Kb(0), k' b {0) = ke{0). Since 7 is a 
spin-2 field that transforms as 7' = je 2tip , this operation is also equivalent to a rotation 
of each ellipticity by 7r/4 with each position vector fixed. 

3'4. Magnification Bias. - In the absence of gravitational lensing, the cumulative 
number counts no(> So) (per solid angle) of background galaxies above the limiting flux 
Sq can be locally approximated around S = So by a power-law form: 

(63) n o( >So) = £^J^ocV* 
with the running power index around S = So 

(64) a^- dl ° gwM>So) >0. 

d\og 10 S 

Note the value of a depends on the luminosity function dn/dL(L, z) of background 
sources (and hence the object type, such as late/early type galaxies and quasars) and 
the observing wavelength. 

Gravitational lensing induces the following conflicting effects known as magnification 
bias (see Refs. [54] and [55]): 

1 Area distortion: Stl 1 = fj,(0)SQ s 

2 Flux Magnification: S -> (j,(6)S. 

The former effect reduces the effective observing area in the source plane, decreasing 
the number of background sources per solid angle; on the other hand, the latter effect 
amplifies the flux of background sources, increasing the number of sources above the 
limiting flux. The net number counts with gravitational lensing is given as 

r°° (ft m 

(65) n(>S )^ dS——=v a - 1 n (>So). 

JSa/LL (J-dUdb 

This implies that, (i) positive bias for a > 1 and (ii) negative bias for a < 1. 
In the weak lensing limit (\k\, \^y a \ <C 1), 

(66) n(> So) w (1 + 2k)"- 1 n (> S ) « {1 + 2(a - 1)k)} n (> S ). 
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The fractional change in the surface number density of background sources is thus given 
as 



Hence the magnification- bias Sn(0) in the weak-lensing limit is a local measure of the 
surface-mass density field, n{0). A combination of the lens distortion and magnification 
can thus be used to break the mass-sheet degeneracy inherent in the shear-based mass 
determination^, 9]. 

4. Cluster Weak Lensing Profiles 

In this section we summarize basic, useful techniques for measuring cluster weak 
lensing profiles, which can be compared quantitatively with various cluster models, and 
can be used to constrain cluster mass and structure parameters. 

4T. Weak Lensing Distortion. - The spin-2 shape distortion of an object due to 
gravitational lensing is described by the complex reduced shear, g = g x + ig 2 = 7/(1 — k), 
which however is coordinate dependent. For a given reference point on the sky, one 
can instead form coordinate-independent quantities, the tangential distortion g + and the 
45° rotated component, from linear combinations of the distortion coefficients g\ and g 2 
as [56] 

(68) g + = -(gi cos2</? + g 2 sm2tp), g x = -(g 2 cos 2ip - g\ sin 2ip), 

where ip is the position angle of an object with respect to the reference position, and 
the uncertainty in the g + and g x measurement is er + = a x = Ogj\f2 = a in terms of 
the rms error a g for the complex distortion measurement. The +-component, g + , is a 
measure of tangential coherence of the shape distortions of background images due weak 
gravitational lensing (see Figure 5 and discussions in §3'3) . On the other hand, the 
x-component, g x , corresponds to divergence-free, curl-type distortion patterns of back- 
ground images. In practice, the reference point is taken to be the cluster center, which 
can be determined from symmetry of the strong lensing pattern, the X-ray centroid po- 
sition, or the brightest cluster galaxy position[13]. To improve the statistical significance 
of the distortion measurement, we calculate the weighted average of the g + } s and its 
weighted error as 



(67) 



Sn = — w — 2(1 — ch)k. 



(69) 



(g+{6 n )) = J , (g x (6 n )) - ' , 

2^ a j 2^ a j 



(70) 



C mn = (Ag + (8 m )Ag+(9 n )) = (Ag x (8 m )Ag x (9 n )) = 8 mn a\{9 m ) 



(72) 



(71) 




j'Gbin n 



where the j runs over all of the objects located within the nth annulus with a median 
radius of 6 n , g+j is the ^-component of the reduced shear estimate for the jth object, 
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statistical weight for the jth object, and C mn is the bin-to-bin covariance error 
matrix of the binned distortion profiles. Here we have used {g a ,i9f),j) — (l/2)(Tg j5 a p5ij 
(a, (3 = 1,2).( 3 ) Several authors adopted the statistical weight u g ,i of the form: u 9i i = 
l/((Tgj + a 2 ), with a being the softening constant variance[57, 58, 40]. In the limit of 
a >• a 2 , this corresponds to the uniform weighting. The case of a = is known as the 

inverse- variance weighting, yielding <j+(9 n ) = <Jx(0 n ) = l/^/Ej(l/^)- ^ Refe. [57] and 
[7] a ~ (a 2 ) w 0.4 was used as the softening parameter for the weight function. 

The tangential reduced shear g+(9) as a function of angular radius is a useful, di- 
rect observable in weak lensing observations, being free from the mass-sheet degen- 
eracy (§3T): It does not require a non-local mass reconstruction (see §3T and 3'3), 
and hence one can easily assess its error propagation as given by equations (70) and 
(71). In particular, the bin-to-bin covariance error matrix C mn is diagonal (assum- 
ing that statistical uncertainties are dominated by random orientations of intrinsic el- 
lipticities). Furthermore, the 45°-rotated g x component can be used as a useful null 
check for systematic effects. To compare observed tangential-distortion profiles with 
cluster mass models, one can use the identify (j + (9)) = k(< 9) — (k(6>))[4][5], where 
(j + (9)) is the azimuthally-averaged tangential component of the gravitational shear 
at radius 8, (k(6)) = (£ m (0))/£ cr it is the azimuthally-averaged k at radius 9, and 
k(< 9) = S m (< 0)/E cr it is the mean lensing convergence interior to radius 9. If the 
projected cluster mass distribution is azimuthally symmetric about its center, then one 
can make an approximation (g + (9)) ~ («(< 9) — (k(9)))/(1 — (k(9))) to compare with 
predictions from azimuthally-symmetric mass models (see Table 1). 

In Figure 9 we show the tangential and 45°-rotated distortion profiles for A1689 at 
Zd = 0.183 (left) and A2142 at — 0.091 (right) derived from Subaru weak lensing 
data[59, 9, 60], along with the respective best-fitting g + (#)-profilcs for the NFW (see 
Ref. [37]) and singular isothermal sphere (SIS) models (see Table 1). In both cases, the 
cluster distortion profiles are better fitted by the NFW model with a continuously steep- 
ening profile. For the nearby, high-mass cluster A2142, the curvature in the distortion 
profile appears highly pronounced, and the SIS model can be strongly rejected by the 
weak lensing data alone. 

4'2. Weak Lensing Depletion. - Lensing magnification fi influences the observed sur- 
face density of background sources, expanding the area of sky, and enhancing the observed 
flux of background sources (§3'4). For red background galaxies at a median rcdshift of 
zj ~ 1, the intrinsic count slope s = —dlog w no(> So)/d\og w So = (5/2)dlog 10 n(< 
mo)/dmo at faint magnitudes too is relatively flat (s < 1), so that a net count deple- 
tion results [54]. Recently the count depletion of red background galaxies due to grav- 
itational lensing has been clearly detected in several massive clusters (A1689, A1703, 
A370, RXJ1347-11, CL0024+1654)[8, 12, 40]. Figure 10 displays the count-depletion 
profile derived from a Subaru £>i? c z'-selected red galaxy sample in the background of 
CL0024+1654 at z s = 0.395 [40]. The error bars include not only the Poisson contribu- 
tion but also the variance due to variations of the counts along the azimuthal direction, 
i.e., contributions from the intrinsic clustering of red galaxies[54, 61] and departure from 



( 3 ) We have ignored the contribution of large scale structure along the line of sight to the 
lensing signal, namely the cosmic shear signal, acting as spatially correlated noise on the cluster 
lensing measurement. See Ref. [35]. 
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Fig. 9. - Azimuthally-averaged radial profiles of the tangential reduced shear g+ [upper panels) 
for the high-mass galaxy clusters A1689 at Zd = 0.183 (left) and A2142 at Zd = 0.091 (right) 
based on deep Subaru weak lensing data. The solid and dashed curves show the best-fitting 
NFW and SIS profiles for each cluster. Shown below is the 45° rotated (x) component, g x . 
Figure taken from Ref. [60] . 



circular symmetry [9] . A strong depiction of the red galaxy counts is shown in the central, 
high-density region of the cluster, and clearly detected out to a few arcminutes from the 
cluster center. The gray-filled region represents the 68.3% confidence bounds for the 
predicted count depletion curve n(9) = n fi s ^ 1 (9) from an NFW model constrained by 
Subaru distortoin data, demonstrating clear consistency between these two independent 
lensing observables. 



4'3. Weak Lensing Dilution. - It is crucial in the cluster weak lensing analysis to 
make a secure selection of background galaxies in order to minimize contamination by 
cluster/foreground galaxies and hence to make an accurate determination of the cluster 
mass, otherwise dilution of the distortion signal results from the inclusion of unlensed 
cluster galaxies, particularly at small radius where the cluster is relatively dense [59] [42]. 

This dilution effect is simply to reduce the strength of the lensing signal when averaged 
over a local ensemble of galaxies, in proportion to the fraction of unlensed cluster and 
foreground galaxies whose orientations are randomly distributed, thus diluting the lensing 
signal relative to the true background level, derived from the uncontaminated background 
population [42]. With a pure red background sample (B) as a reference, one can quantify 
the degree of dilution for a galaxy sample (G) containing JVcl cluster galaxies and Nbg 
background galaxies in terms of the strengths of the averaged tangential shear signal 
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Fig. 10. - Number-count profile of Subaru Bi? c z'-selected red galaxies (squares) in the back- 
ground of C10024+1654 at Zd = 0.395. The triangles show the counts without the mask cor- 
rection due to cluster members and bright foreground objects. A slight dip at 9 = 2' — 3' 
in the depletion profile corresponds to the contribution of the known substructure (see Figure 
6), which is also seen in Subaru distortion data. The gray-filled region represents the 68.3% 
confidence bounds for the predicted count depletion curve from an NFW model constrained by 
our Subaru distortion analysis, demonstrating clear consistency between these two independent 
lensing observables. Figure taken from Ref. [40]. 



(g+(9)) as [42] 

N bg + N C l _ (9 { + B \0)) (Dds/Ds)[% d 



(73) 1 + 8 d {6) 



Nbg (gf\6)) (D ds /D S )[ B) , 



where (Dd s / ' D a ) Za ^ Zd , s are averaged distance ratios for respective background popu- 
lations^); if the two samples contain the same background population, then S d = 
(.9+ ) / (.9+ ) ~ !■ The degree of dilution thus varies depending on the radius from 
the cluster center, increasing towards the cluster center. Medezinski et al. [42] found for 
their green galaxy sample ([V — i']~E°/so-o 3) containing the cluster sequence galaxies in 
A1689 that the fraction of cluster membership, A*cl/(^Vbg + Nql), tends ~ 100% within 



( 4 ) Background samples can be generalized to include foreground field galaxies with z s < Zd 
and D ds /D s (z s ) = 
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R < 300 kpc^ 1 . 

Broadhurst, Takada, Umctsu et al. (see Ref. [59]) proposed to use a sample of red 
background galaxies whose colors are redder due to large fc-corrections than the color- 
magnitude relation, or red sequence, of cluster member galaxies. These red background 
galaxies are largely composed of early to mid- type galaxies at moderate redshifts[42]. 
Cluster member galaxies are not expected to extend to these colors in any significant 
numbers because the intrinsically reddest class of cluster galaxies, i.e. E/SO galaxies, are 
defined by the red sequence and lie blueward of chosen sample limit, so that even large 
photometric errors will not carry them into such a red sample. This can be demonstrated 
readily, as shown in Figure 11, where we plot the mean tangential shear strength (g + ) 
of A1689, averaged over a wide radial range of 1' < 9 < 18', as a function of color limit 
by changing the lower color limit progressively blueward. Here we do not apply area 
weighting to enhance the effect of dilution in the central region. Figure 11 shows a sharp 
drop in the lensing signal at A(V — i') ^ 0.1, when the cluster red sequence starts to 
contribute significantly, thereby reducing the mean lensing signal. At A(V-i') ^ 0.1, the 
mean lensing signal of the red background stays fairly constant, (g + ) ~ 0.143, ensuring 
that our weak lensing measurements are not sensitive to this particular choice of the 
color limit. 

Recently, this empirical background selection method has been generalized by Medczin- 
ski, Broadhurst, Umetsu et al. [62] to incorporate and combine all color and positional 
information in a color-color (CC) diagram. This CC-selection method has been success- 
fully applied to Subaru imaging observations of several massive clusters[62, 40]. 

On the other hand, the dilution of the lensing signal caused by cluster members can 
be used to derive the proportion of galaxies statistically belonging to the cluster, or 
the cluster member fraction, by comparing the undiluted background distortion signal 
with the radial distortion profile of color-magnitude (or color-color) space occupied by 
the cluster members, but including inevitable background galaxies falling in the same 
space [42]. This technique allows the light profile of the cluster to be determined in a way 
which is independent of the number density fluctuations in the background population, 
which otherwise limit the calculation of the cluster light profiles and luminosity functions 
from counts of member galaxies. The resulting light profile can be compared with the 
mass profile to examine the radial behavior of M/L[42, 62]. 

4'4. Aperture Mass Densitometry. - The observed image distortion of background 
galaxies can be directly used to derive the projected gravitational mass of clusters. The 
aperture mass estimate within the angular radius 8, M^(< 9), in terms of the tangential 
component 7+ of the gravitational shear can be expressed as 



using the aperture-mass densitometry, or the so-called "^-statistic" [65] , defined as 



where 9 out (> 0) is the (constant) outer background radius, (j+(9)) — k(< 9) — (k(9)) is 
an azimuthal average of the tangential component of the gravitational shear at radius 9, 
and k = S m /S cr it is the mean convergence. Errors on £(9) are calculated by propagating 
the rms errors a + (9 n ) (equation [71]) for the tangential shear measurement. Equations 



(74) 



M c (< 8) 



(75) 
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Fig. 11. - Top panel: mean distortion strength averaged over a wide radial range of 1' < 6 < 18' 
of A1689, done separately for the blue and red galaxy samples. No area weighting is used here 
to enhance the effect of dilution in the central cluster region. Shown are the measurements of 
the tangential component (g+) with open squares, and those of the 45deg-rotated component 
(<?x). Bottom panel: the respective numbers of galaxies as a function of color limit, contained 
in the range 1' < 6 < 18' in the red (right) and the blue (left) samples. This figure is taken from 
Ref. [9]. 



(74) and (75) show that the cluster mass can be measured from the galaxies ellipticity 
within the annulus bounded by 9 and 9 ou t located just outside the mass to be measured. 
In the weak lensing regime where \k\, \j\ <C 1, (-f+(9)) is observable: (g + (6)) ps (j + (9)). 

As revealed in equation. (75), the £— statistic yields a mean convergence interior to 
9, subtracted by the mean background within the annulus between 9 and 6* out , 7c(< 
9) — k(6, 9 ont ). Hence, as long as k(9, 9 out ) <C k(< 9), the enclosed mass within 9 can be 
obtained by multiplying £ by the area 7rS cr it(Z?d^) 2 - The inner radius 9 can almost be 
arbitrarily chosen to obtain the aperture mass interior to 9, as long as the weak lensing 
approximation is valid. Obviously, the aperture mass is smaller than the enclosed mass by 
a negative compensating mass that serves to remove the contribution from a background 
uniform mass sheet, and the degree of deviation depends on how steep the density profile 
is. 

The ^-statistic is a circular aperture mass estimator. It gives rise to some errors for 
irregular and merging clusters. Nevertheless, this effect has been estimated less than 
10% [63]. The choice of the parameter 9 out may also affect the cluster mass estimate in 
a practical application. For example, a small 9 out will generate large Poisson noise since 
the galaxy number for a shear estimate within the annulus bound by 9 and 9 out is not 
sufficiently large. On the other hand, if 9 ou t is too large, the cluster mass measurement 
can be contaminated by neighboring clusters in projection. For a projected lens system, 
(7+) is produced not only by the cluster itself but also by the projected neighboring 
clusters and large scale structure. 

Clowe et al. [64] proposed a variant of aperture mass densitometry [65], by introducing 
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two parameters to specify the annular background region, of the form: 

Cc(0) = 2 / dln6'( 1+ (6')) + / dln0'< 7+ (0')) 

(76) =K(<0)-K(0 inn A ut ), 

where #i nn and 9 out ($out > $inn > 0) are the inner and outer radii of the annular 
background region in which the mean background contribution, Rb = K(9- lnn , 9 out ) , is 
defined. This cumulative mass estimator is often referred to as the Cc-statistic. This 
cumulative mass estimator subtracts from the mean convergence R{9) a constant Rb for 
all apertures 9 in the measurements, thus removing any DC component in the control 
region 9 = [9[ nn , 9 out ]. We note that the Rb is a non-observable free parameter. The 
^-statistic will be particularly useful for wide-field imaging observations, such as with 
the Suprime-Cam on the Subaru telescope (34' x 27') and the Megacam on the CFHT 
(1° x 1°), in which one can identify a background region well outside of the cluster region. 

4'5. Weak Lensing Convergence. - Umetsu & Broadhurst [9] have developed a non- 
parametric method for reconstructing the one-dimensional K-profile utilizing the ( c - 
statistic measurement. Unlike strong- lensing based boundary conditions [8], this method 
utilizes an outer boundary condition on the mean background density Rb to derive a 
K-profile (see Schneider & Seitz 1995 for an alternative method for a direct inversion of 
the mass profile). ( 5 ) For a given boundary condition Rb, the average convergence k(< 6) 
is estimated as 7c(< 9) = ( c (9) + Rb- Then, we define a discretized estimator for k as 

(77) K{6 m ) - <C(0m+l) - <C(M + Kb, 

where 

n, m - 1 (°m\ 2 m _ 1 ( 8m+l \ 2 

and 9 m is the weighted center of the mth radial bin bounded by (9 m , 6 m +i) (see Refs. [9] 
and [40]). The error covariance matrix C mn of K m is expressed as 

C mn = a?a%Ci +hn+1 + - ar«2C£,„ +1 - «c£ +1 ,„, 

where C£„ = (5£ m 5( n ) is the bin-to-bin error covariance matrix of the aperture densit- 
ometry measurements which is calculated by propagating the rms errors a+{9 m ) for the 
tangential shear measurement (see §4T). In the non-linear regime, however, the 7+(0) is 
not a direct observable. Therefore, non-linear corrections need to be taken into account 
in the mass reconstruction process. In the subcritical regime (i.e., outside the critical 
curves), the 7+(0) can be expressed in terms of the the averaged tangential reduced 
shear as (g + (9)} w -f + (8)/[l — n{9)\ assuming a quasi circular symmetry in the projected 
mass distribution. This non-linear equation (76) for Cc(#) can be solved by an iterative 
procedure, which is outlined in Ref. [40]. 



( ) The mass-sheet degeneracy is inevitable in any mass reconstruction method based solely on 
the shearing effect. 
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In the left panel of Figure 12 we show cluster surface mass density profiles £ m (0) 
for A1689[9] and CL0024+1654[40] as reconstructed from combined Subaru weak-lensing 
(r ^ 200 kpc) and HST/ACS strong- lensing (r 200 kpc) observations. The joint mass 
profiles for the clusters continuously steepen out to their virial radii, and are well fitted 
with single NFW profiles, which provide a good description of the equilibrium den- 
sity profiles of collisionless DM halos in cosmological iV-body simulations. The right 
panel of Figure 12 shows two-dimensional marginalized constraints on the NFW model 
parameters (e v ir,M v i r ) derived for CL0024+1654 from the joint mass profile shown in 
the left panel. This figure demonstrates that combining strong and weak-lensing in- 
formation (grey contours) significantly reduces the uncertainties on the profile param- 
eters. Such non-parametric mass profiles are also useful when comparing the total 
matter distribution with cluster properties obtained from other wavelengths and/or 
approaches [66, 67, 68, 69]. 




Fig. 12. - Left: Cluster surface mass density profiles for A1689 (green triangles) and 
CL0024+1654 (red squares) reconstructed from strong lensing (HST/ACS: r ^ 200 kpc) and 
weak lensing measurements (Subaru: r ^ 200 kpc). The curvature of these profiles demon- 
strates that both strong and weak lensing data are required to make an accurate measurement 
of the cluster mass structure parameters, such as the halo virial mass, M v j r , and concentration, 
Cvir- Also shown as solid curves are the best-fit NFW models, with a continuously steepening 
density profile, for A1689 (black) and CL0024+1654 (gray). Right: Joint constraints on the NFW 
model parameters (c v i r , M v i r ) derived for CL0024+1654. The red and blue contours show the 
68.3%, 95.4%, and 99.7% confidence limits for the inner strong lensing and outer weak lensing 
data, respectively. Combining strong and weak-lensing information (gray contours) significantly 
reduces the uncertainties on the profile parameters. Figures taken and modified from Umetsu 
& Broadhurst [9] and Umetsu, Medezinski, Broadhurst et al. [40]. 



Appendix A. 

Multiple Lens Equation 

We may discretize the cosmological lens equation (19) by dividing the radial integral 
between the source (x = Xs) and the observer (x = 0) into A-comoving boxes (N — 1 
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lens-planes) separated by a comoving distance of A%. The angular position {n > of a 
light ray at the nth plane (n < N) is then given by 

(A.i) - = - V Kxn - x P ) v±v 

where ip p is the effective lensing potential of the pth lens plane (p — 1,2,..,N — 1): 

(A.2) V> P =-r/ rfx*- 

The Jacobian matrix is expressed as 

(A.3) A<"> :=_■ = /- £ Xn )jyW AW = I — 

^ P =i 

where ff^g = d 2 tp p /dx ± ' a dx ± '^ {a,fi — 1,2) is the Hessian matrix, and g{xpiXn) is 
the effective lensing distance for the pth lens plane. In general (without the Boron 
approximation), the 2x2 Jacobian matrix can be decomposed into the following form: 

(A. 4) A a fi = 5 al } - ip a0 = (1 - n)8 a[i - 7i cr 3 - 7 2 cti - ioJcr 2 

where k = (i/>n + i>22)/2, 7i = (V'n ~ ^2i)/2, 72 = (^12 + ip2i)/2, and w is the net 
rotation w = (V'12 — V'2i)/2. The Born approximation A*- 1 '- 1 = / leads to the symmetric 
Jacobian matrix. 

* * * 
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